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Abstract 

We propose that, in M = 4 Super Yang-Mills theory, correlation functions 
of certain 1/8 BPS Wilson loops and local operators inserted on a S 2 in space- 
time may be computed in terms of analogous observables in the "zero-instanton" 
sector of 2d Yang-Mills theory. The Wilson loops are mapped to the standard 
Wilson loops of the 2d theory, as recently conjectured, while the local operators 
are mapped to powers of the 2d field strength. We give several perturbative 
checks of the correspondence, and derive from 2d Yang- Mills a two-matrix model 
for the correlator of a local operator and a Wilson loop of arbitrary shape. We 
show that the strong coupling planar limit of the two-matrix model precisely 
agrees with a string theory calculation in AdS§ x S 5 . 
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1 Introduction 

Supersymmetric Wilson loops in the four- dimensional M = 4 super Yang-Mills theory 
have been extensively studied over the past years. One specific motivation is that, 
in certain cases, they may provide examples of physical observables which are non- 
trivial and yet exactly calculable. In particular, one may obtain this way interesting 
quantitative tests of the duality to type IIB string theory in the AdS§ x S 5 background 

In 0, 0] an exact result has been conjectured for the circular maximally supersym- 
metric 1/2 BPS Wilson loop operator: its expectation value can be computed using 
a Gaussian Hermitian matrix model. This conjecture has passed many subsequent 
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tests, in particular it agrees with all the available calculations in the dual string the- 
oryB In [16] the Gaussian matrix model has been derived using localization of the 
four-dimensional path integral to supersymmetric configurations. 



In 17H19I] a large class of interesting, generically 1/16 BPS, Wilson loops has been 
found. Those loops live on a three-sphere S 3 in Euclidean space-time R 4 . Impos- 
ing some restrictions on these Wilson loops on S 3 one gets various 1/8, 1/4 and 1/2 
BPS Wilson loops. In particular, restricting to the equator two-sphere S 2 C S 3 , one 



gets generically 1/8 BPS Wilson loops. In [17H19l| it has been conjectured that the 
expectation values of 1/8 BPS Wilson loops on S 2 are exactly captured by a purely 
perturbative calculation in the two-dimensional bosonic Yang-Mills theory on S 2 . In 
two dimensions, the preferred gauge choice is the light-cone gauge, since then there are 
no interactions. The conjecture of 17H19I] then implies that the 4d expectation values 



gaugj 



should be equal to the sum of the ladder diagrams of 2d Yang- Mills in light-cone 
This prescription is not equivalent to the exact 2d bosonic Yang-Mills theory 
but instead to a "truncation by hands" of all non-zero instantons on S 2 24l-l26l|. In 
271 . |28| the conjecture has been supported at A 2 order for a single Wilson loop on 



S 2 (in the case of the "two-longitudes" loop [27J, |28| as well as "wavy-latitude" loops 
(j^l ). where A = gyM^ * s the ^ Hooft coupling constant for Yang-Mills theory with 
SU(N) gauge group. In particular, in 27J, |28( it was found that even when there are 
non-trivial M = 4 SYM interacting Feynman diagrams, the final result agrees with the 
ladder diagram computation in 2d YM in light-cone gauge. However, for a connected 
correlator of two latitudes on S 2 , [28( found a discrepancy at order A 3 between the Feyn- 
man diagrams in M = 4 SYM and the ladder diagrams of 2d YM in light-cone gauge. 
Subsequently, in [29j and [30( several new tests in support of the original conjecture 
have appeared. In particular it was shown in 29j that invariance under are preserving 
diffeomorphisms holds at strong coupling, as implied by the conjecture, and a Gaussian 
two-matrix model for the connected correlator of two Wilson loops on S 2 was derived 
from the exact solution of 2d YM 29J 30j. Its strong coupling limit agrees with the fact 
that there are no connected supersymmetric string worldsheets joining the two loops 
29j |. and the first subleading corrections to the saddle point at strong coupling have 



been shown to agree in 30| with the exchange of light supergravity modes between 
the two worldsheets. Further, in 30| a Feynman diagram calculation in M = 4 SYM 
at order A 3 , in the limit of one shrinking loop, was shown to be consistent with the 
two-matrix model. 



In 31J the localization framework was used again to understand the relation be- 



2 See [fjj for a review of earlier work, and [7- 15 1 for a partial sample of more recent relevant work. 
3 More precisely, this is an Euclidean version of the light-cone gauge, defined by A s — 0, where A 
is the 2d gauge field and z, z are complex coordinates on S 2 . 



2 



tween the four-dimensional M = 4 SYM and the two-dimensional theory on S 2 , where 
the interesting Wilson loops live. Using localization, one gets naturally a Lagrangian 
formulation of the conjectured 2d theory, which turns out to be closely related to 2d 
Hit chin/ Higgs- Yang-Mills theory 32j-|34|, and also a natural explanation of the pre- 
scription to truncate the 2d instantons in the 2d YM conjecture. The localization 
computation in j3l| for 1/8 BPS loops was not completed at the same level of rigour 



as m 



161 ] for 1/2 BPS loops. One still needs to evaluate the one-loop determinant 



for the field fluctuations in the directions normal to the localization locus. However, 
there are many reasons to believe that such determinant is trivial in the M = 4 theory, 
and then the localization 3l| would support the original conjecture 17H19I] . Moreover, 
the localization framework allows to establish a complete correspondence between all 
observables in the M = 4 SYM which share a number of certain superconformal sym- 
metries and observables of the 2d theory. In particular, one immediate consequence 
is establishing the 2d description of certain local operators on S 2 which share some 
common superconformal symmetries with the relevant 1/8 BPS Wilson loops. These 
operators are chiral primaries equipped with an explicit space-time dependence, of the 
form tr(— ^-H^b)" 7 , where x % are coordinates on M 3 on which the two-sphere x l x l = r 2 
is embedded, $j are the three scalars which couple to the 1/8 BPS Wilson loops, and 
$s is any of the remaining three scalars. Note that the definition involves an identifica- 
tion of a 5*0(3) subgroup of the i?-symmetry group with the SO (3) rotating the x l, s, as 
it is natural in the construction of the 1/8 BPS Wilson loops of jl7H19|. Actually, these 
local operators on S 2 are a special case of a more general class of protected operators 
on IR 4 which was recently studied in 35(, where the preserved supersymmetries and the 
non-renormalization properties of their correlation functions were investigated. When 
an arbitrary number of such local operators is inserted on S 2 , the system preserves 
4 superconformal supercharges. On the other hand, the combined system of Wilson 
loops and local operators on S 2 preserves 2 common supercharges, which is sufficient 
for the localization of |3l| to be applicable. In particular, it follows that these local 
operators are mapped in the 2d theory to insertions of powers of the YM field-strength 
(or more precisely its Hodge dual). 

In this note we make several detailed weak and strong coupling tests of this cor- 
respondence involving the chiral primaries on S 2 , and our results support the 2d YM 
conjecture. In particular we study, to leading order in perturbation theory, the correla- 
tor of a local operator and a Wilson loop and obtain agreement with the corresponding 
computation in 2d YM. Further, from summing up the ladder diagrams of 2d YM in 
light-cone gauge we derive a two-matrix model for the the exact correlator of a local 
operator and a Wilson loop. This can also be written as a complex matrix model, 
and then one can see that our results imply as a special case the original conjecture 
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of 36] for the exact correlator of a 1/2 BPS circular loop and a chiral primarjQ. We 
solve the two-matrix model in the planar limit, and show that at strong coupling it 
exactly matches the corresponding string theory calculation in AdS§ x S 5 , using the 
explicitly known string solutions for the 1/4 BPS latitude and the 1/4 BPS 

two- longitudes [nl loops. 



Hence, until now, all available results essentially support the original 2d YM con- 
jecture, except the result in 28] on the connected correlator of two Wilson loops at 



A 3 order. In an attempt to explain such discrepancy, in [30] doubts were raised as 
to whether the light-cone gauge prescription A z = is equivalent to the Hermitian 
two- matrix model 29|, [30| which captures the zero-instanton sector of 2d YM on S 2 . In 
appendix we show that the ladder Feynman diagrams in the light-cone gauge A s = 
for the correlator of two latitude Wilson loops on S 2 are actually in precise agreement, 
to all orders, with the Feynman diagrams of the Hermitian two-matrix model, so the 
discrepancy in 28] still remains unsolved. 

The paper is organized as follows. In Section [2] we set up our notations and con- 
ventions, we explain the 2d description of the local operators implied by localization, 
and we study the supersymmetries preserved by local and Wilson operators on S 2 . In 
Section [3] we give our perturbative checks of the proposed 2d-4d correspondence. In 
Section H] we derive the two-matrix model from 2d YM in light-cone gauge and solve it 
in the planar limit. In Section [5] we present the string theory calculation of the corre- 
lator between a local operator and a Wilson loop. Finally, in the Appendix we collect 
some notes about 2d YM in light-cone gauge and the equivalence with the Gaussian 
matrix models derived from the zero-instanton sector. 



2 Preliminaries 

2.1 Notations and conventions 

The M = 4 SYM action on R 4 with the standard flat metric is 

Ssym = — ^- / d 4 x [\tiF^F^ + tvD^ A D^ A + . . . j , (2.1) 
9ym J \ z J 

where [i = 1, . . . , 4 are space-time indices and A = 1, ... ,6 are SO(6)r indices. Here 
we use conventions such that the covariant derivative is D = d + A, the curvature 
is F^ u = [D^jD,,], and all fields take value in the Lie algebra of the gauge group, 

This was extended to the 1/4 BPS circular loop in [371 ] . See also [(5, 3^- 42 1 for additional related 
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work 
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= A^Ta, $ = $ a T a , e.g. in the anti-Hermitian matrices for the U(N) gauge group. 
The anti-Hermitian generators satisfy tr T a Tb = —\5 a b- Hence the action may be also 
written as 



S 



SYM 



' f d 4 x ( l -FlFl + DMDM + ...]. 



(2.2) 



The 1/8 BPS Wilson loops of jl7l4l9|| are located on a sphere S 2 of radius r defined 

= 1,2,11 



as £4 = 0, Yli=i x l = r2 1 an d they couple to three of the six scalars, $ 



W R {C) 



-f- tr R Pexp (b (Aj + ie ijk ®i—)dx 
d R Jc r 



x* 



(2.3) 



where dn denotes the dimension of the representation R. According to the conjecture 
of |17H19||. this supersymmetric operator is mapped to the standard Wilson loop on the 
same contour C and representation R in the two-dimensional Yang-Mills theory with 
action 



S 



SYM 



1 r #*vg QfyFr 



where 



2 

9id 



did 



2nr' 



(2.4) 
(2.5) 



We will use the notation A, F for the fields of the two-dimensional theory. In 
particular 31] we have A = A + i * 2 d ®t where * 2 d is the Hodge stai@ on 5* 2 and 
$ t is the two-component one-form obtained from the components of the 4d field $j 
"tangent" to the S 2 , see eq. ( I2~3l) . 



2.2 Localization for local operators on S' 



In 3l|] it is shown that the 4d M = 4 SYM path-integral localizes to a 2d theory 
on S 2 , namely to the constrained Hitchin/Higgs- Yang-Mills theory, or conjecturally 
to the zero-instanton sector of the standard bosonic Yang-Mills, which we denote as 
aYM theory (here "aYM" stands for "almost Yang-Mills", in view of the fact that 
contributions of the unstable instantons are dropped). 



The localization computation [31[ implies that certain local observables inserted on 
the same two-sphere where the 1/8-BPS Wilson loops are located, are also mapped to 

J Hcre we use the conventions in H3E3. These differ from the conventions used in by a relative 
sign in the scalar couplings. 

6 Our conventions for the Hodge star are such that in flat space with metric ds 2 — dx\ + dx\ we 
have *2ddx 1 — dx 2 , *2ddx 2 — —dx 1 , and the orientation on S 2 is the standard orientation of flat space 
when we use the stereographic coordinates. 
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the two-dimensional theory. We briefly explain this fact in the following, and refer the 
reader to 31 for more details on the localization calculation. 



Choose one of the three remaining scalars which do not couple to the Wilson loops 
on S 2 and denote it $5. In 31] it is shown that at the localization locus one has 
$ B = 0, whikfl 



dT^t 



1 



(2.6) 



where <3> t is the "tangential" one-form obtained from $j as explained above, and $ n = 
Yli=i denotes the component "normal" to the S 2 . Here d*2 d = *2ddA*2d and 
dA = d + A. Next, in 31] we find another localization equation 



d A *2d ®t = iF, 



(2.7) 



which relates the 4d field $i to the 2d field F. Hence, combining (I2.6p . (12.71) and 
$b = 0, we see that at the localization locus we have the relation 



+ i$ B = ir * 2d F. 



(2. 



The field $ n + i§B on S 2 is Q-closed, where Q is the fermionic symmetry used in the 
localization computation (this is one of the two superconformal supersymmetries shared 
by the Wilson loops and local operators on S 2 , see Section [2731) . Hence localization is 
applicable to operators which are gauge invariant functionals of this field. Therefore, 
the results of [3JJ imply that the operator Oj(x) = tr($ n + i^bY for x G S 2 in the 
M = 4 SYM theory is mapped to the operator tr(ir *2d F) J in the 2d aYM, and hence 
the correlation function of any number of such operators and any number of Wilson 
loops are mapped to the corresponding correlation functions in the two-dimensional 
aYM theory. 

In this paper we aim to explicitly compute some correlation functions of this type 
and compare with the strong coupling limit using the dual AdS$ x S 5 description. 



2.3 Supersymmetry 

We now show that the supersymmetric Wilson loops on S 2 in (12. 3p and the local 
operators Oj(x) = tr($„ + i^b) J share two preserved supercharges. 

It is convenient to use the notation of M = 1 SYM in lOd, and split the 10 Dirac 
matrices as Tm = (7 m ,Pa), where /x = ,4, A = ,6, 7^ are space-time 



7 Literally in 3l| one finds the equation d*2 d §t — —&n- We ch ang ed the sign because of differ- 
ent conventions in the definition of the Wilson loop (|2.3[) versus [31(, and also added the explicit 
dependence on r. 



6 



gamma-matrices and pA are the S0(6)r Dirac matrices. The combined variation of 
the bosonic fields under Poincare and superconformal supercharges can be written as 



SAp = , S^a = ippAe , (2.9) 

where ip is the gaugino and 

6 = 60 + ^7^1, (2.10) 

where Sq, t\ are 16-component spinors corresponding respectively to the Poincare and 
superconformal supercharges. 

Let us first review the supersymmetries preserved by the Wilson loops (12.31) . To 
simplify notations, we will set r = 1 throughout this section. The loops live on x A = 0, 
x l x l = 1, so we split the 4d gamma- matrices as (74,74). Moreover since they only 
couple to three of the scalars, we write pa = (pi, p±, P5, pe), where the index i = 1, 2, 3 
is identified with the space-time 3d vector index, and p±, p$, p§ are rotated by 5*0(3) b C 
SO(6)r. The supersymmetry variation of the loop (12.31) then yields 

8W oc x 3 x k (jjk^i + i£ijkPi^o) ~ x l jix J x k (7 jfc eo + i£ijkPiti) ■ (2.11) 

Therefore the variation vanishes for arbitrary loops provided that 

7 ifc ei + ie ijk pie = . (2.12) 

One can eliminate for example ei from these equations to obtain the following condi- 
tions 

(7y+p i i)eo = 0, i,j = 1,2,3. (2.13) 

These are three consistent equations, but only two are independent since the commu- 
tator of any two equations gives the remaining one. With two independent projectors, 
we are left with 4 independent components of eo- Using any of the equations in (12.121) . 
the superconformal spinor e\ is completely determined in terms of eo as 

ei = -^Pi23eo- (2.14) 

So the conclusion is that the loops (12.31) preserve 4 combinations of Poincare and 
superconformal supercharges. 

We now turn to the local operators on S 2 , and let us choose = $4 to be concrete. 
In this case, the supersymmetry variation yields (see also 351 ]) 

50j(x) oc x l (pi€ + iptliti) ~ x l x 3 lj {Piti + ipilito) . (2.15) 

The variation vanishes independently from the insertion point x % if 

Pi^o + ipiliti = . (2.16) 
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As before, we can proceed by eliminating e± from these equations, which yields the 
constraints 

(7y + Pii)eo = 0, i,j = 1,2,3, (2.17) 

which are exactly the same conditions found for the Wilson loops above. So again we 
have 4 independent solutions for e . The superconformal spinor can be now determined 
in terms of 6q from any of the (I2.16p . and the result is 

ei = -«7iPiP4e • (2.18) 

Therefore the conclusion is that the local operators on S 2 preserve 4 supercharges, but 
these are not the same 4 supercharges preserved by the Wilson loops. 

To see whether the local operators and the loops share some supercharges, we 
should impose that (12.141) and ( 12. 18ft are simultaneously satisfied. This yields the 
further condition on e 

(p 23 + 7iP 4 ) eo = . (2.19) 

One can now see that this condition is consistent with the three equations fl 2 . 1 3 j) . as 
commutators of any two of the four equations ( 12. 13ft - (12. 19ft either vanish or produce an 
equation in the same set. Therefore there are three independent projectors and hence 
2 independent solutions for eo- The spinor t\ is given in terms of eo using either (I2.14p 
or ( 12.18p . so we conclude that the combined system of any number of local operators 
and any number of Wilson loops on S 2 preserves 2 supercharges. 



3 Explicit perturbative checks 

In this section we give some explicit perturbative checks of the correspondence Oj(x) <-> 
tr(ir *2d F) J between 4d and 2d theory. 

One could always use conformal invariance to take the S 2 to have unit radius r — 1, 
but we have chosen not to do so to keep dimensions of the relevant quantities in the 
2d and 4d theory more transparent. For simple book-keeping let us summarize the 
dimensions of the relevant quantities in terms of unit of length [L\ 

[x] = [r] = L 1 [g 2 YM } = L° [$''] = IT 1 [F^] = L~ 2 [W(C)] = L° [Oj\ = L~ J 
[A=L° [gl d ] = L- 2 [F- ZZ ] = L° [* 2d F}=L~ 2 , 

(3.1) 

where z is the complex coordinate on S 2 , see below. 

In our conventions the 4d propagators for the gauge field in Feynman gauge and 
for the scalars are 
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In the 2d theory, we work with complex coordinates on S , with metric given by 

2 Ar 2 dzdz 

ds = — — -g , (3.3) 
(1 + zz) 

where r is the radius of the sphere. This is related to the standard metric in polar 
coordinates by z = tan |e^. For 2d perturbative calculations, it is convenient to use 
the "Euclidean light-cone" gauge defined b}|§ 

A- z = 0. (3.4) 

In this gauge there are no interactions and the 2d YM action becomes simply 

Sym 2 = "TtV / d 2 z^g z %A a z g zz d- z A a z . (3.5) 

We use notations d 2 z = dz A dz = 2idx A dy for z = x + iy, and ^fg = —ig ZZ i so that 
d 2 Zy/g is the conventional volume form on S 2 normalized as 



d 2 z ^g = 47rr 2 . (3.6) 
The gauge field propagator is [18 § (see also the Appendb 



i /./,,,, \ ulr ,„i, 1 1 z-w 



A a z (z)A z (w) ) = ^—5 a \ — - — - . (3.7) 

zK ' zK 'I IX 1+zzl+wwz-w 

This satisfies^ 

(V9(9 Z *) 2 dz (A a z (z)A z (w))) = 6 ab 6 2 (z - w) . (3.8) 

92d v x 1 ' 

On S 2 there is no ambiguity in the propagator for the kinetic term ( 13. 5ft . A quick 
explanation is that the kinetic term operator is the square of the Dolbeault operator d, 
which maps (1, 0)-forms on S 2 ~ CP 1 , which are sections of the bundle 0(— 2), to the 
(1, l)-forms. However, the bundle 0(—2) does not have holomorphic sections, hence 
there are no zero modes, and the propagator is well defined. 



8 The components of the 2d gauge field A z and A s are treated as independent. This choice of gauge 
is consistent for perturbative calculations, but it cannot capture the non-perturbative corrections since 
there are no classical solutions (instantons) that satisfy this gauge. 

9 In [l8| the gauge field propagator is given in a certain generalized Feynman gauge, in which both 
A z and A z are propagating, and the A z propagator is related to the one in the "light-cone" gauge 
used here by a factor of 2. 

10 Our convention is such that / d 2 z5 2 (z — w)f(z, z) — f(w, w) for any f(w, w). 
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We can also explicitly show that there is no ambiguity by solving the equation 



d- z (p- 2 d- z G z (z,z))=0 (3.9) 

on C and analyzing the behaviour at infinity Here we denoted p = (1 + zz)~ l . From 
lpT5]) we get 

d- z G z (z,z) = p 2 f(z), (3.10) 
where f(z) is an arbitrary holomorphic function on C. Solving (13.101) . we get 

G z (z,z) = -^—f(z)+g(z), (3.11) 
1 + zz 

where g(z) is again an arbitrary holomorphic function on C. Now we must require that 
G z at z — > oo decreases at least as fast as z~ 2 in order for the solution to be smooth 
on CP 1 . Indeed, if we make a coordinate transformation to the coordinate z = 1/zso 
that the point z = oo maps to the point z = 0, we get G z dz = G z dz = —G z j§. Asking 
G z to be finite at z — 0, we see that G z decreases at least as z~ 2 at z —>■ oo. However, 
the solution (13.111) implies that G z (z, z) decreases not faster than z -1 , unless both f(z) 
and g(z) vanish. Therefore, the equation (I3.9P does not have non-zero smooth solutions 
on CP 1 . The actual solution (13.71) has the correct asymptotics z~ 2 at z — ► oo, hence it 
exists and is well defined globally on CP 1 . 

For later convenience, we also write down the explicit expression for the scalar dual 
to the 2d field strength in these coordinates 

*2d F = ^-F- zz = ^- (1 + zz) 2 d- z A z . (3.12) 
3.1 Correlators of local operators 

In the 4d theory, the tree level 2-point function of the elementary fields making up 
the local operators Oj(x) on S 2 is (here and in the following we pick $ B = $ 4 to be 
concrete) 



(($„ + ^ 4 ) a (a;)($ n + z$ 4 ) 6 (y)) = -xy(^(x)^(y)) - (nW&S)) 

9YM_fiab ( X ' V/ r2 ~ ^ _ 9YM jafe 



An 2 \ (x — y) 2 J Sn 2 r 



(3.13) 



where we have used x 2 = y 2 = r 2 since the operators are inserted on S 2 . Thus we 
see that correlation functions between the local operators OAx) = tr($„ + z$4) J are 



position independent at tree level. This was also observed in 35[ (for the more general 
operators on M 4 of which our operators on S 2 are a special case), and it was argued that 
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position independence holds true even at the quantum level by using a Ward identity 



which follows from the preserved supersymmetries. It was further argued in |35j, based 
on the results of ji^ . that correlation functions involving an arbitrary number of Oj(x) 
do not receive quantum corrections (up to possible instanton corrections) and thus can 
be obtained by simply doing Wick contractions with the free propagator (I3.13p . The 
position independence implies that the n-point function can be computed by a Gaussian 
multi-matrix model. It is not difficult to see that the following n-matrix model correctly 
captures the gauge group combinatorics 



J nMM 



[dXf\ ■ ■ ■ [dX«] %) ^ 1 ' ' , (3.14) 



where Xk = X^T a are matrices in the Lie algebra of the gauge group. Due to position 
independence, the tree-level n-point functions are then given by 

(0 Jl (x 1 )Oj 2 (x 2 )---0 Jn (x n )) = (tTX(HrX^..-tTX^) nMM , (3.15) 

where the correlation function on the right-hand side is taken in the n-matrix model 
(I3.14p . Inverting the quadratic form in the matrix model action, one can see that there 
are only propagators between different matrices, as appropriate for normal-ordered 
operators Oj(x), and the normalization of the action in (13.141) is chosen to match the 
propagator (13.131) . 

Let us now turn to the 2d theory. To check the correspondence between Oj(x) and 
ti(ir *2d F) J , let us compute the 2-point function 

(ir * 2d F a (z) %r * 2d F b (w)} = ^(1 + zzf(\ + wwfl,d- z A^z)d^~A h z (w)) . (3.16) 
Using the light-cone gauge propagator (13. 1\ we get 

(ir * 2d F a (z) ir * 2d F b (w)) = 5 ab (°k- l -9k(l + zz) 2 5 2 (z - w^j . (3.17) 

The 5-function piece does not matter if we consider correlation functions of normal 
ordered operators inserted at distinct points. If we ignore the 5-function term, we 
then see that this propagator agrees with (I3.13P provided we identify the 2d and 4d 



2 

couplings as g\ d = — f^?, as implied by the conjecture of [171, ll8( for Wilson loops on 
S 2 and by the localization calculation of 3l|. Since in the light-cone gauge the 2d YM 



action is free, correlation functions of ti(i* 2d F(z)) J inserted at arbitrary points will be 
also given, perturbatively, by just doing contractions with the free propagator (13.171) . 
The combinatorics for the contraction of gauge indices is the same as the one for the 
4d local operators, so we can conclude that 

{Oj x { Xl ) ■ ■■0 Jn {x n )) Ad = (tr [ir * 2d F(zS) ' ■ ■ -tr (ir * 2d F{z n )^j Jn )° 2 j Bt , (3.18) 
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where the equality holds for normal-ordered operators inserted at distinct points. The 
non-renormalization arguments of f35| and the localization argument we propose in 
this paper suggest that the relation (13. 18[) should hold to all orders in perturbation 
theory (possible instanton contributions on the 4d side are not in principle excluded). 

3.2 Correlator of a local operator and a Wilson loop 

We now consider the correlation function of a local operator and a Wilson loop on S* 2 , 
to leading order in perturbation theory. We will start by considering the correlation 
function involving a single elementary field <3> n + 2$4, and let us suppress the gauge 
indices for the moment. In the 4d theory, we wish to compute 

(^$,(x)+i$ 4 (^ W{C)). (3.19) 

Expanding the path-ordered exponential in the Wilson loop, to leading order we have 
to evaluate 

£(A j + ie ijk ^W) =i^£dyi^0^ 2 . (3.20) 

It is not difficult to compute this integral for arbitrary loop, see also [18]. Define 6 to 
be the angle between x l and y % . If we denote by d<p the one-form orthogonal to dd, 
then we get 

/ dyj ^V k = _I # sin2g = _ / #C0S 2 1 . (3.21) 
Jc (%-y) 2 Jc V 2(l-cos0) Jc 2 K J 

For r = 1 this is equal to half the area of the region of S 2 enclosed by the loop and not 
containing the point x, up to a choice of +/— sign which depends on the orientation of 
the loop relative to the point x. If we denote by S + ,S~ the two regions of S 2 singled 
out by the loop, see Figured], such that S + is the one containing the north pole, and if 
we take the loop to run counterclockwise with respect to the north pole, then we can 
summarize the result as 

/ r l \ r „ k ( j 9ym A 2 rz C+ 

((^W + ^Wj/^ + fe***^ - | +< |^; i£S _ • (3-22) 

where A\ and A 2 are the areas of S + and S~ respectively. It can be seen that this is 
precisely the behavior expected in the 2d Yang-Mills theory, in particular we see that 
the result only depends on the area and not on the shape of the loop. Moreover, it is 
almost independent from the position of the local operator, it only depends on whether 
the operator is inserted "inside" or "outside" the loop. 
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Figure 1: The correlator of a Wilson loop W(C) and a local operator 
Oj(x) on S 2 . The curve C divides the S 2 into two regions which we 
denote as S + and S~. Perturbatively, we find that the correlator only 
depends on whether the operator is inserted in the S + or S~ region. 



Let us now carry out the analogous calculation in the 2d theory. Here, to first order 
in the coupling, we have to compute 

(ir * 2d F(z) j dwA w ) . (3.23) 

Since we are working at first order in perturbation theory, we can essentially assume 
that we are in the abelian theory. Then by Stokes' theorem we can write 



A= F , (3.24) 

c Js+ 

where S + is the region enclosed by the loop and containing the origin of the complex 
plane (i.e. the north pole). Using (13.121) we get 

(ir *2d F(z) <p dwA w ) = - [ d 2 w 2r (ir * 2 d F(z) ir * 2d F(w)) 
Jc Js+ I 1 + wwy 

9l<f [ d 2 w ; 2 ( 1 _ 2m (i + W uj) 2 5 2 (z-w)) 



47T Jq+ (1 + WW) 



-i§±A u zeS- 



(3.25) 

which indeed agrees with the 4d result (I3.22j) by virtue of the identification g 2 d = 
~9ym /2vrr s As a double-check of factors of i and signs, let us also do this computation 
in a different way using Cauchy's theorem. Using the invariance under area preserving 
diffeomorphisms of 2d YM, we can for simplicity consider a latitude on S 2 , which is a 
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circle of radius tan-^ in the z, z coordinates (here 6$ is the latitude angle). Then 

(ir * 2d F(z) <f dwA w ) = + zzf I dw (d- z A z {z)A z {w)) 

JC lr J\w\=tan^ 



cos — (b dw h tan - 

2vr 2 ,/ H=tan £a \w-z 2 w(w - z) 



-47rcos 2 f , |z|<tanf 
;4£ . /i^oiT.a «a i„i ^ + OT1 eo 



-^•47rsin z f, |z| > tan f, 

(3.26) 

where we have used that w = r^/w for a circular contour of radius ro, and in the last 
equality we have used Cauchy's theorem. As expected this agrees with f )3.25p . since 
A 2 = Airr 2 cos 2 y and A\ = Anr 2 sin 2 for a latitude circle. 

It is straightforward to extend the above calculations to the case of the correlator 
(Oj(x)W(C)) 4d (or equivalently ( tv(ir * 2 d F) J W(C) ) ) for arbitrary J. To first order 

\ /2a! 

in perturbation theory, we consider a Feynman diagram obtained by joining the local 
operator to the Wilson loop with J "local-to-loop" propagators as computed in eq. 
(13.20P (or eq. (13.251) ). Since the product of the J-propagators is symmetric in the 
exchange of the vertices, path ordering of the exponential in the Wilson loop does not 
matter. Considering only the planar contractions we then get 



1 ( N\ J ( .g\ M I J d Sjj h x l y \ , 2 ./+i 



(3.27) 



(O^WiC))^-^-) ^iSULfV^ + o i9YM 



N(j-i)\\(=£ry , xe s- + [ h 

where A = Qym^ ls ^ e Hooft coupling. In the first line the factor 1/N comes from 
the normalization of the Wilson loop (we take R to be the fundamental representation), 
the factor (—N/2) J comes from the contractions of gauge group generators, and the 
factor of J counts the number of planar diagrams obtained from cyclic permutations. 
At the next orders in perturbation theory, one decorates the above Feynman diagram 
by ladder "loop-to-loop" propagators, as well as internal interaction vertices and loops 
in the case of the 4d theory. In the 2d YM theory in light-cone gauge, on the other 
hand, there are no interactions and the full perturbative result is given by summing 
up ladder diagrams. If the conjectured 4d-2d correspondence is correct, the sum of the 
light-cone ladder diagrams in the 2d theory should be equal to the sum of all Feynman 
diagrams (including interacting diagrams) in the 4d SYM theory. In the next section 
we will write down a Gaussian two-matrix model which computes the sum of the 2d 
light-cone ladder diagrams. 
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4 A matrix model for the correlator of Oj and Wr{C) 



Let us start by considering for simplicity the case of a latitude on S 2 . The result can 
then be generalized to arbitrary loops by using the invariance under area preserving 
diffeomorphisms of the 2d YM theory. In the case of the latitude, it is especially simple 
to sum up ladder diagrams, because the "loop-to- loop" propagator is a constant. Let 
us parameterize the loop as z(r) = r e lT , where r = tan y and 8 is the latitude angle. 
Then the propagator between two points Ti,r 2 along the loop is, using (13.71) . 

(z 1 A M (z 1 )z 2 A M (z 2 )) = S - (1 + rg)2 =^ (27r)2 A . 

where in the last equality we have written the result in terms of the areas singled out 
by the loop, and A = A\ + A 2 . Since the loop-to-loop propagator is a constant, path 
ordering is not important and the calculation of the previous section directly shows 
that the correlator is independent of the insertion point of the local operator. Hence, 
specializing to the case in which the local operator sits in S + , we can effectively consider 
the insertion point to be the north pole, i.e. z = 0. Then the local-to-loop propagator 
is also a constant, see eq. (13.261) 

(ir , M f (0) zA\(z)) = W^ffi^ = • (4.2) 

This allows us to easily sum up the 2d light-cone ladder diagrams in terms of the 
following two-matrix model, which we write directly in terms of the 4d coupling using 
the relation g\ d = -\g\ M 

t r y 2 --^-tr xy) 



Z 2MM = / [dX\[dY\ e ^ > , (4.3) 



where X and Y are N x N Hermitian matrices. The propagators in this matrix model 



arc 



which match respectively the loop-to-loop and the local-to-loop propagators, eq. (14. ip 
and (14. 2 p (to make the comparison, use g\ d = —\g YM and recall that in our conventions 
the U(N) generators satisfy (T a )*(T a )| c = — \5\8^, and that we integrate over the loop 
t G [0, 27r]). The (YY) propagator vanishes as required by normal ordering of the 
local operator. Therefore the correlators of interest are given by the two-matrix model 
correlation functions 

(tr(ir *2 d F(z)) J W R (C)) ~ inSt = (tiY J -^-ti R e x ) , (4.5) 

\ / 2d \ d R / 2MM 
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where the label "O-inst" is to remind us that the sum of the 2d light-cone ladder dia- 
grams does not capture the non-perturbative contributions of the unstable instantons 
on S 2 . In the above we have assumed that the local operator is inserted in the region 
including the north pole, and we will do the same in what follows. If the operator is 
inserted in the complementary region, one simply exchanges the roles of A\ and A 2 
and includes an additional factor of (—1) J due to orientation. 



The arguments of 17J, ll8|], the localization calculation of 311] and the identification 
tr(ir *2d F) J «-> Oj(x) lead us then to propose that the same two-matrix model may 
exactly capture the correlators in the 4d M = 4 SYM theory 



(Oj(x)W R (C)) 



1 



-ui 



J 2MM 



[dX][dY] trY J 



1 



d 



tr R e x e 



2g YM 



-try 2 



-tiXY) 



R 



(4.6) 

Note that for general loops this matrix model is not equal to the sum of the ladder 
diagrams in the 4d theory, since the combined gauge-scalar ladder propagator is not 
equal to the 2d light-cone propagator, but rather to its real part [l7J, [l8|. If the loop-to- 
loop propagator is a constant, which happens in the case of the 1/4 BPS latitude, then 
the matrix model is indeed equal to the sum of 4d ladder diagrams, and the proposed 
relation to 2d YM would imply that the sum of interacting diagrams should vanish. 
It may be possible to explicitly check this^to leading order in perturbation theory, by 
adapting to our case the results of 0, 36-38, 43]. For more general loops, on the other 
hand, one should check whether the sum of ladder and interacting diagrams in 4d is 
equal to the 2d light-cone ladder diagrams. In the case of the 1/4 BPS loop made 
up of two arcs of longitudes (l7| . this has been successfully checked in 27, 28] for the 
case of the expectation value of the Wilson loop. It should be possible to extend those 
calculations to the case of the correlation function with a local operator. 

As a consistency check of the two-matrix model, notice that if we do not insert the 
local operator, then we can integrate out Y exactly and we end up with 



(W R (C)) 



1 



■id 



J MM 



1 



[dX]^-tr R e x e 2Al ^" 
dR 



i — 

YM 



trX 2 



(4.7) 



which is precisely the matrix model proposed in 17|, ll8j to exactly capture the expec- 
tation value of the 1/8 BPS Wilson loops on S 2 . In the case of a loop at the equator, 
A\ = Ao = A/2, it reduces to the well-known matrix model for the 1/2 BPS circular 
loop 

In the next subsection we will solve the two-matrix model in the planar limit. 
Remarkably, the strong coupling behavior precisely agrees with the dual string theory 
calculation we present in Section [51 
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4.1 Large N solution 

It is straightforward to solve the two-matrix model (14. 6 ft in the large N limit, following 
45j (see also 29J for more details on the application of the formalism of 45] to the 



Gaussian two-matrix model). Since the action is Gaussian, it is also straightforward to 
solve the matrix model at finite N, but we do not do this here. For simplicity, in the 
following we will restrict to the case of a Wilson loop in the fundamental representation. 

For a general Gaussian two-matrix model with action 

W/ = Ar(^trX 2 + |trr 2 - C trXr) , (4.8) 
the planar two-point resolvent 



w(*i,*2) = (tr-_ ^tr — - ) . (4.9) 



ci,>mi 



is given by 

u(zi, z 2 ) = -d Zl d Z2 log (1 - yi{z 1 )y 2 {z 2 )) , (4.10) 
where yi(zi) and y 2 {z 2 ) are two resolvent functions 

Vi(zi) = ^ ( Zl + yj z i~ 4 ««i^) > 2/2(^2) = 7^ (z2 + \J~4~ 4aa 2 ^ • (4.11) 
Here the parameters a, a\ and a 2 are related to the parameters in the action by 

2 C °2 a l iq\ 

a = -, ofi = — a, a 2 = — a. (4.12) 

ai«2 — c c 

For the two-matrix model (14. 3p . we get 

2 , A 2 kAiA 2 . . 

a =zA-prr, Aaai = — -5 — A, a 2 = . (4-13) 

As it is clear from the definition (14. 9p . to extract the correlators (tr Y J tr e x ) we should 
expand the resolvent in inverse powers of z 2 



uizi, z 2 



J2^n(zi)z 2 n ~ 1 , (4.14) 



n=l 



extract the term with n = J, and perform the inverse Laplace transform on the Z\ 
variable 

Wj(h) = (f dziuj^e* 1 ' 1 . (4.15) 



2ttz . 
Then 

<trF J tre* lX ) = Wj{t{) . (4.16) 
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By expanding (I4.10p . we get 

u J (z 1 ) = Ja J yi J - 1 d zl y 1 . (4.17) 
Using the explicit expression for yi(zi), the inverse Laplace transform yields 

Wj{t x ) = Ja J 2 /j(2t lv ^r) , (4.18) 

where Ij(x) is a modified Bessel function of the first kind 

Ij{x) = — [ d<f)e iJtl> e xcos * . (4.19) 
2vr Jo 

Putting everything together, our final result for the correlator in the planar limit is 
then 
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.,1. x \ i J J\" 2 r' f A 2 \' , I [AMM 



<°^ W ^ = ( bY N tie )-N^-{fj MV^ A| ' (420) 
Expanding this expression at weak coupling, we get 

(OAx)W(C)) = V (l + + 0{ y) I , (4.21) 



iV(J-l)! V^ 2 / V ^ 2 J + ! 
which agrees, as expected, with the leading order perturbative result ( 13. 27ft . 

Let us now extract the strong coupling asymptotics of (14.201) . To compare to the 
string theory calculation, it is somewhat more convenient to work with the normalized 
local operators, see Section [5] 

dj(x)=2 J /\^y^=Oj(x). (4.22) 

At strong coupling it is also natural to normalize the correlator by the Wilson loop 
expectation value, which is given by jil. 17. 18 



(W(C)) = -^h(^) , A' = A -^\ . (4.23) 



The normalized correlator then reads 



oMmo) _ (4 , 4) 



(W(C)} N \AJ \aj j^v^ 7 ; 
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Using the asymptotic expansion^*] 



IM 1 _^ pl+ J t --»f + 3 + ... ; (425) 



2a; 8a; 2 
we get the following large A expansion of the correlator 



(W(C)) N \AJ \AJ V V A 

(4.26) 

The leading term is in precise agreement with the dual string theory on AdS§ x S 5 , as 
shown in Section [5j 



4.2 Equivalence with the complex matrix model 

The two-matrix model (14.61) can be written in complex notations, so that one can 
see the equivalence with the complex matrix model which was first proposed in 36 
to compute the correlator of a circular 1/2 BPS Wilson loop and a chiral primary 
operator, see for related work and [37[ for the extension to the circular 1/4 

BPS Wilson loop. 

First we rewrite the matrix model (14.61) as 

, 1 Y ^! tr (zl±L Y {2X+^Y)) 



(Oj(x)W R (C)) M = [ [dX][dY] trY J ^ti R e x e 

"2MM J UR 



(4.27) 



Next we introduce new complex matrices 



~ lAl -Y Z = 2X + ^Y (4.28) 



and change variables in the matrix integral (I4.27P from X and Y to Z and Z@ The 
Jacobian of this change of variables is a number independent of X and Y, which only 
changes the overall normalization of the partition function. Since we normalize the 
correlation functions in the matrix model by the partition function, this Jacobian does 
not matter. Therefore we get 

— / [dZ dZ) tr Z J — tr R e^ z+z ^ e~^^ . 

(4.29) 

11 Here we only include the contribution of the dominant saddle point at x — » oo. 

12 Note that Z and Z defined in (|4.28[) are not complex conjugate to each other, but the number 
of degrees of freedom is still the same. Perturbatively, we can always rotate formally the integration 
contour so that Z and Z are complex conjugates. 
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This agrees with the complex matrix model of 36| and its extension to the 1/4 BPS 



circular loop 37J. The % in the factor (ir/Ai) J takes into account that $ in our 
conventions is valued in the Lie algebra and thus is an anti-Hermitian matrix for the 
SU(N) theory, and the r/A\ comes from the finite-distance propagator between our 
local operator and the Wilson loop. The coefficient in the action of the complex matrix 
model in the case of a latitude at polar angle 6 reduces to the familiar [371 



A 2 



2g 2 M A x A 2 



ti ZZ 



9YM sin 



tr ZZ. 



(4.30) 



Note that in earlier works one often considers the case in which the local operator 
is inserted at large distance from the Wilson loop, so that the correlator is used to 
extract the coefficient of the chiral primary in the operator product expansion of the 
Wilson loop [46]. However, suppose that we insert the local operator at the north pole 
of S 2 . Then we can make a conformal transformation which maps the north pole to 



infinity, and our results on the matrix mode 
have shown that the 2d YM- conjecture of 



for Oj and W(C) will still hold. Hence we 
17H19| on the Wilson loops on S 2 , refined 



here to treat the local operators (see also |3JJ), implies as a special case the complex 
matrix models suggested in [36|, |37|, |39( for the OPE of circular 1/2 (1/4) BPS Wilson 
loops and chiral primaries. 



5 Strong coupling computations 

The single trace chiral primary operators of M = 4 SYM take the form 

C AlA2 ... Aj ti^ A ^ M ---^ , (5.1) 

where A%, . . . , Aj = 1, . . . , 6 are SO(6)r indices, and the tensor Ca x a%-Aj is symmetric 
and traceless in each pair of indices. These operators preserve half of the Poincare 
supersymmetries, and they have protected conformal dimension A = J. In particular, 
an operator of the form 

tr(u-$) J , (5.2) 

with u a complex six- vector satisfying u 2 = 0, is a chiral primary operator. Note that 
the operators of interest in this paper, Oj = tr ($„ + £$4) J , take this form and are 
therefore chiral primaries. However notice that the complex vector u is taken to be 
spacetime dependent u = (x l /r, i, 0, 0) in this case. 

The AdSs x S* 5 dual description of the chiral primary operators is well established 
from the early days of the AdS /CFT correspondence. A chiral primary operator of 
dimension J corresponds in the bulk to the Jth Kaluza-Klein mode on S 5 of a certain 
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supergravity scalar field which is a linear combination of the fluctuations of the metric 
and RR 4-form potential along the S 5 directions^. This was first worked out in 50], 
based on the analysis of 51], and we refer the reader to the original literature for more 
details. Let us denote the relevant lOd scalar field by s(x, y), where x^ = (z, x), x e 1R 4 
are coordinates on AdS$ with Poincare metric 



ds 2 = — (dx 2 + dz 2 ) 



(5.3) 



and y are coordinates on S 5 . Then we can expand the field s(x,y) in S 5 spherical 
harmonics 

s(x,y) = J2s J (x)Y J (y), (5.4) 



where Y J (y) are scalar spherical harmonics satisfying 



V 2 F J 



-J(J + 4)F J . 



(5.5) 



Note that spherical harmonics on S 5 transform in the same representation of SO(6)r as 
the chiral primaries (15.11) . If we parameterize the S 5 in terms of flat 1R 6 coordinates Q A 
satisfying G 2 = 1, then the spherical harmonic corresponding to (15.11) and satisfying 
(I5.5P is just 

Y J (Q) = C AlM ... Aj Q M Q M ■■■Q A \ (5.6) 

or Y J = (w • 6) J for the operator in (15.21) . The supergravity equations of motion imply 
that the 5d scalar fields obey [50| [51 



V 2 s J = J(J-A)s J . 



(5.7) 



By the standard AdS /CFT dictionary, a scalar field in AdS^, with m 2 = J(J — 4) 
is dual to an operator of conformal dimension J, which is identified with the chiral 
primary operator at the boundary. In the spectrum one finds only modes with J > 2, 
consistently with the fact that the dual gauge theory has gauge group SU(N) (chiral 
primaries with J = 1 do not exist due to tracelessness of the SU(N) generators). 

The equation of motion for the 5d field s J (x) with a source located at the boundary 
can be solved in terms of the bulk-to-boundary propagator 



s J (z, x) 



d 4 x'G j(z, x; x')sq(x') 



(5. 



13 This description is appropriate as long as J <C N. Chiral primaries of large conformal dimension, 
J ~ N, are dual to spherical D3 brane probes wrapping an S 3 inside AdS§ or S 5 . These are the giant 
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gravitons of 
backreacted geometries of 49 1 



For dimension J ~ N the most appropriate dual description is given by the 
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where Sq is the source at the boundary, and the propagator is given by 

z J 

Gj{z,x-x')=Nj f, (5.9) 

V ; (z 2 + (x - x') 2 ) J 

where Mj is a normalization factor. According to the standard AdS /CFT procedure, 
correlation functions of the chiral primary operators are obtained by plugging the above 
solution into the supergravity action and functionally differentiating with respect to 
the source. It is convenient to chose the normalization constant Aj so that the chiral 
primaries corresponding to Sq have unit normalized two-point functions. This depends 
on the absolute normalization of the type IIB supergravity action, see [50(] , and fixes 



the normalization constant to be 46 



Ny/J 



The correspondingly normalized operators on the gauge theory side take the for 

Oj(x) = 2 J / 2 H/^= tr (u ■ $) 3 , (5.11) 

which satisfy 

(Oj^O*,^)) = 6jj ' . (5.12) 
{xi - x 2 ) 

Here we have assumed that u is a constant six- vector satisfying u 2 = and u ■ u* — 1. 
In the following string calculation we will adopt the same normalization conventions 

for our operators Oj{x) = tr yy&i + i^A j on S 2 with x-dependent w-vector. This 
means that we consider the normalized operators 

dj(x) = tr (j^ + t* 4 ) ' • (5.13) 

Their two-point function satisfies (note that here we do not want to take the complex 
conjugate) 



(dj(x 1 )dj / (x 2 )) = JJ ' 2J • u(x 2 )) 



- x 2 ) 

j {x\ ■ x 2 /r 2 -l) J = (— i) 5jj,. 



(xi - x 2 

Sjj' / , i ( — 1 



j (5-14) 



(2(r 2 -xi-x 2 )) j v ' V2r 2 

From the supergravity point of view, the factor (u(xi) ■ u{x 2 )) J arises from the inte- 
gration of the two spherical harmonics (u(x{) ■ 0) J and (u(x 2 ) ■ 0) J ' over S 5 , while the 
denominator comes from the bulk-to-boundary propagator in AdS 5. 



14 The perhaps unfamiliar factor of (— i) J is due to our conventions that the fields $^4 are anti- 
Hermitian. 
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5.1 Correlators of local operators and Wilson loops 



In the string theory dual, the correlator between a local operator Oj(xq) and a Wilson 
loop is obtained by computing the amplitude for the process in which the supergravity 
mode dual to Oj(x ) is emitted from the insertion point x at the boundary and then 
absorbed at a point on the string worldsheet dual to the Wilson loop [46J (see also 
0] for a review), see Figure [2J This requires to work out how the scalar field s(x,y) 
couples to the string worldsheet. The Polyakov action for the AdS$ x S 5 string reads, 
in conformal gauge 

S = J d 2 a (g^d a x^d a x u + g mn d a y m d a y n ) , //, v = 1, . . . , 5 , m, n = 1, . . . , 5 , 

(5.15) 

where g^, g mn are the AdS 5 and S 5 metrics respectively. To linear order, the supergrav- 
ity mode s(x,y) is related to the metric fluctuations h^, h mn around the background 
g^u, g m n as follows [40|, |46j, [50 

6 4 
V = —^Js(x,y)g^ + — -V(jid v) s(x,y), 

hmn 2,7 s(x, y)gmn ■ 

In the first line V^c^) stands for the symmetric traceless part 

V {u d u) s(x, y) = V u d u s(x, y) - -g^V^x, y) . (5.17) 

To linear order, the coupling of the s(x,y) field to the worldsheet is simply obtained 
by inserting the first order fluctuations (15.161) into the Polyakov action 

S {1) = ^,\ d2 ° (K v d a x»d a x» + h mn d a y m d a y n ) . (5.18) 

The field s(x, y) can be written in terms of the boundary source as explained in the 
previous section 

s(x,y) = J d A x'Gj{z,x-x')si{x')Y J {y) , (5.19) 

and the normalized correlator is then obtained by functionally differentiating with 
respect to the source 

(Oj{x )W{C)) _ SS (i) 

(w(c)) -~64W)- (5 - 20) 

The calculation now proceeds analogously to [371 . |46j. One important difference, 
however, is that in those works the local operator was inserted at large distance from 
the Wilson loop, which yield a considerable simplification in the form of the bulk-to- 
boundary propagator. In our case, we wish to insert the local operators on the same 
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Figure 2: The correlator (Oj(xo)W(C)) at strong coupling. The plane 
represents the S 2 G R 4 at the AdS$ boundary, on which the Wilson loop 
C resides. The operator Oj(x ) is inserted at a point x G S 2 , in either 
of the regions S + or S~ . The wavy line represents the supergravity mode 
dual to the local operator, which propagates from the insertion point Xq 
to a point on the string worldsheet dual to the Wilson loop. 



S 2 on which the Wilson loops live, so we must work exactly in the separation between 
the local and Wilson operators. 

For the explicit calculation, it will be convenient to work with the flat embedding 
coordinates Q A on S 5 satisfying 2 = 1. The spherical harmonic corresponding to 
Oj(xo) then can be written as 

Y J (Q) = (u(x ) ■ 0) J = + iQ 4 ) , z = l,2,3. (5.21) 

The string solutions dual to the 1/8 BPS Wilson loops on S 2 reside on a AdS-s x S 2 
subspace of AdS 5 x S 5 , where the AdS 3 G AdS 5 is defined by x % x % + z 2 = r 2 , = 0, 
and the S 2 G S 5 is parameterized by the unit 3- vector 0*, while 4,5,6 = 0, see (l7| . 
Hence the 4 in (I5.2ip can be effectively dropped in the following. 

After carrying out the functional differentiation with respect to the source, we get 

^T^T = IT-, l d2 ° f -2 J( f ~ 1} Gj(x»; x )\d a x»d a x>* + -J—d a d a Gj(x»; x )) Y 
osq(xo) Am a' j V J+l z l J+l / 



'o 1 

6S$} 



—2J [ d 2 ad a B A d a e A Gj(x fl ;x )Y J (e). 
not 1 J 



(5.22) 



8sq(xq) Ana' 

To obtain the second term in the first line, we have used the fact that 

daX^d^V^Gj = d a d a Gj , (5.23) 

which can be proven by noting that the string worldsheet satisfies the equations of 
motion 

d a d a af + T^ p d a x u d a x p = . (5.24) 
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To further simplify (I5.22p . it is convenient to integrate by parts the second term in the 
first line. This produces a boundary term which can be neglected since Gj goes to zero 
sufficiently fast at the boundary (at least if the insertion point Xq does not touch the 
Wilson loop, which we assume to be the case), plus a term proportional to 

d a d a Y\Q) = - i f ' d { *® X Q )2 daQ - Jd a e A d a e^j y j (&) . (5.25) 

Here we have used the equations of motion for the Q A 

d a d a Q A + d a Q B d a Q B Q A = 0. (5.26) 

Then the AdS§ and S 5 contributions in (I5.22p can be combined to give our final 
result 

& - ^ («• + m) - >H Gjyj • 

(5.27) 

where we have used the relation a' = 1/ y/X. By plugging in an explicit string solution 
and evaluating the above integral one obtains the strong coupling prediction for the 
correlator (15.201) . 



Let us start by considering the string solution dual to the 1/4 BPS latitude [171. 143 
The corresponding Wilson loop on S 2 is depicted in Figure [HH- The explicit solution 
is given bjo 

i tanhcxocosr 2 tanh <To sin r s 1 4 

x = r , x = r , x = r , x = 0, 

cosh a cosh a cosh <7o 

z = r tanh <7o tanh a , (5.28) 

^1 cost ~ sinr „ , , , 

= r^' = H n. = tanh a + a) . 

cosh(o" + a) cosh(o" + erj 

Here the range of the coordinates is < a < oo, < r < 2tt, and the parameter 
<t is related to the latitude angle on S 2 by tanhcr = sin^Q. Inserting this solution 
into (I5.27I) . together with the explicit expression for the propagator given in ( 15.9ft . 
one obtains a somewhat complicated expression which we do not explicitly report 
here. Remarkably, we have verified by direct numerical integration that the result 
does not depend on the precise position x G S 2 of the local operator, but only on 
whether the operator sits "inside" or "outside" the loop. This is precisely the behavior 
expected from the gauge theory analysis and the relation to 2d YM. The integrand 



15 Here we consider only the stable solution. One can carry out the analogous calculation for the 
unstable solution found in 4c 
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Figure 3: The 1/4 BPS latitude (a) and two-longitudes (b) Wilson loops. 
The corresponding dual string solutions are given respectively in eq. 15.281 
and eq. 15.321 



simplifies considerably if we take the insertion point to be the north or south pole: 
xq = (0, 0, ±r, 0). By doing a change of variables £ = tanher we obtain 



_ .- J(±1)J 2- j / 2 - 1 (j - i)^^ 1 r 1 dc (e + tio\ J e (i + a + ^ 



(5.29) 

where £o = tanhao = sin^o, and the choice of sign corresponds respectively to the 
insertion point at the north and south pole. Notice that the £ integral is convergent 
since J > 2, and the integration is elementary 



1 , r ( e+&o \ J toO-+z 2 )+% i ( e+tto 

So our final result is 



J 



(5.30) 



Oj(x )W(C) J r -j ( (cos |) J+i( sin } XoeS - 

2 v JX \ 



(W(C)} N [(-l) J (sin|)- 7+1 (cos|)-- /+1 , x eS- 

(5.31) 

We see that this is in precise agreement with the matrix model result (14.261) . since 
At /A = sin 2 f and A 2 /A = cos 2 f . 

Another string solution which is explicitly known is the one corresponding to the 
1/4 BPS loop made out of two half longitudes The corresponding Wilson loop on 
S 2 is depicted in Figure [3b. In conformal gauge, the string solution is (l7| [2~* 

i a sin aa sin a + cos aa cos a 2 a cos aa sin a — sin aa cos a 
x = r , , x = r- 



cosh vl — a 2 r cosh \/l — a 2 r 

3 



— r tanh Vl — a 2 T , x 4 = , z = r Q (5.32) 

cosh vl — a 2 r 

6 1 = sin aa , 2 = cos aa 6 3 = , 
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where the range of the coordinates is— oo<r<oo,0<a<7r, and the parameter a 
is related to the opening angle between the two longitudes by 5 = ir(l — a). Plugging 
the explicit solution into (I5.27P we again find that, rather non-trivially, the correlator 
is (almost) independent of the insertion point xq as expected, and the final result is 

(dASo)W(C)) _ r -J f (1 + qffll - a) =4P , x E S + 

(W(C)) - N V \(-l)^(l- a )^(l + a)^, x e S~ ' 

(5.33) 

This is again in agreement with the matrix model result, since for this loop one has 
At/ A = (1 - a) /2 and A 2 /A = (1 + a) /2. 
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A Notes on light-cone gauge vs Gaussian matrix 
models in 2d YM 

In this appendix we collect some notes about 2d Yang-Mills theory on S 2 in the "Eu- 
clidean light-cone gauge" A s = 0. In particular, we first re-derive in an alternative 
way the gauge field propagator by relating it to the two-point function of the field 
strength. Further, in the next section we directly prove the equivalence between light- 
cone gauge Feynman diagrams and the two-matrix model for the connected correlator 



of two latitudes derived in 29] 30] from the zero instanton sector of the 2d YM theory. 

All formulas below are only about the 2d theory, so in this appendix we do not 
have to distinguish the 2d fields from the 4d fields and we do not write tilde for the 
2d fields. The 2d coupling constant is denoted by g 2c i- For simplicity, we will also take 
the S 2 to have unit radius. 

We use complex coordinates z, z on S 2 . The metric for radius r = 1 takes the form 

, 9 Adzdz , . , 

ds = m A -i 

1 + zzY 



so we have 



t. ■ , 2 , g zz = 0, fe = 0. (A.2) 
1 + zz Y 
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The volume form on S 2 is 

H = -idz A dzg zz = , (A.3) 

(1 + \z\ 2 ) z 

where z = X\ + 1x2 and d 2 x = dx\ A dxi = — \dz Adz = — \d 2 z. 

The 2d YM action in the gauge A z = (we skip the Lie algebra indices assuming 
contractions where needed) 



S = -^ ^F- ZZ F ZZ (A.4) 

explicitly takes the form 

S = / -idz A dzFl{2p 2 )- 1 = [ d 2 xp- 2 F 2 z , (A.5) 
where we have introduced 

' = T+li' (A ' 6) 

We now wish to represent the correlation functions of A z in terms of correlation 
functions of F Sz = d z A z . 

We can change variables in the path integral from A z to d z A z . The Jacobian of 
this change of variables is trivial, but in the integration domain over F zz we need to 
explicitly project out the zero modes F zz dz A dz = c/i where c is a constant, because 
such modes are not in the image of d : f2 1,0 — > fl 1 * 1 . Hence, from the free action (1A.5I) 
we immediately get the correlation function of the free fields F zz 

(F- zz (x)F- zz (x')) = -g 2 2d (p(x) 2 6 2 (x - x') - r p^ x)a P(*) 2 f>V) 2 ) > ( A - 7 ) 



which, of course, agrees^j with (13.171) since 

d 2 xp{x) 2 = 7T. (A.8) 



Next we express A z in terms of F zz . Using the fact that 

d- z - = tx5 2 {x) (A.9) 

z 

one easily gets 

Mz) = -L / d 2 u^-. (A.10) 
2m J z — u 



16 In our conventions J d 2 zS 2 (z) — 1, and J" d 2 a:(5 2 (a;) = 1, and d 2 z — dz A dz — 2id 2 x hence 
5 2 (z) = ^S 2 (x) for z = Xi + ix2- 
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The relation flA.101) makes sense on CP 1 if J d 2 uF Uu = 0. 

Now, using (lA.lOp and OA .70 we represent the propagator for A z by means of aux- 
iliary integrals over it-planes 



(A z {z)A w {w)) 



(2m)' 



d u d u 



2„./ 



z — uw — u 



7 (F Uu (u)F m (u')) (A.ll) 



The ^-function term in ( 1A.7I) removes one integral and for the remaining integration 
we use 

d 2 u 



1 



1 



1 



2i (z -u)(w — u) (1 + |u| 2 ) s 



-7T- 



z — w w 

+7T- 



(1 + |z| 2 )(l + \w\ 2 ) Z-W (1 + \w\ 2 ) (1 + \z\ 2 ) 

(A.12) 

This identity can be shown by doing the integral over circles \u\ — const using residues 
and then integrating over \u\. The contribution of the second term in ( 1A.7I) is obtained 



using the integral 



— d 2 u 
2m 



p(u) 



~r~~T\2- (A.13) 

Z — U 1 + \Zf 

Hence, the contribution of the second term in flA.7[) is precisely cancelled by the second 
term in (1A.12j) and we get 



(A z (z)A w (w)) 
which, of course, agrees with ( 13. 7(1 . 



92d 



z — w 



7T 1 + \z\ 2 1 + \w\ 2 Z — W ' 



(A.14) 



A.l Connected correlator of two circular Wilson loops on S 2 

Here we explicitly compute the Feynman diagrams for the connected correlator of two 
latitude Wilson loops on S* 2 using the propagator (IA.14j) for the gauge fields. We prove 



directly the equivalence of the light-cone gauge with the Hermitian two-matrix model 



of 291 30 



We consider two concentric contours, the first contour C\ given by \z\ — r\ and the 
second contour C 2 given by \z\ = r 2 . We assume that r\ < r 2 . 

We denote points on C\ as Wi and points on C 2 as Zi. There are three types of 
propagators (here we use the relations Wi = rf/wi and Z\ = r\jzi): 

From C\ to C\ contour: 



9id ~ r i dwi dwj 

7T (1 + r 2 ) 2 Wi Wj 



(A^dwiA^w^dwj) = ; 1 — -— ± . (A.15) 
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From C*2 to C2 contour: 

(A z (zi)dziA z (zj)dzj) = 9 f^y d f d -^ ■ (A.16) 

From C2 to C\ contour: 

(A z (z)dzA w (w)dw) = ik ~* ] ( z-rl/r\ W \ ^ = 

7r (1 + r{) (1 + r^) V wz(z — w) J 

—r\ 1 I 1 r\jr\ — 1 >^ / w \ n 

(1 +rf) (1 + r|) I wz wz ~[^ z ' 

Now consider a ladder Feynman diagram where we take n points on the contour 
C2 and n' points on the contour C\. A typical diagram is depicted in Figure HI There 
are two type of points 2%. The points z k of the first type are connected by propagators 
to points u>/(fc) on the contour Ci, where I(k) labels the point which connects with z k . 
The points of the second type on C2 are pairwise connected with each other. We denote 
the set of the first type as T(2, 1) and the set of the second type as a disjoint union 
of T(2,2) and J(T(2,2)). In other words, the set T(2, 2) contains a half of the points 
of the second type, and the set J(T(2,2)) contains the remaining half. Let T 2 i be the 
number of points in T(2, 1), and T22 be the number of points in T(2, 2). Analogously, 
the T(l, 1) denotes the points on C\ connected to the points in J(T(1, 1)) on Ci, and 
the connection map is denoted by J(k), i.e. we say that a point w k for k e T(l, 1) 
connects to a point wj(k)- Clearly, n + nl = 2(T 2 2 + T21 + Tn). 

As usual, these Feynman diagrams arise from expanding each Wilson loop in the 
correlator (Wr{C2)Wri{Ci)) in powers of the gauge field 

00 

W R (C) = l + J2w { n\C), (A.18) 

n=l 

and taking all Wick contractions with the gauge field propagator. For example the 
term with n points on C2 is 

W { x\C 2 ) = d R l dz n --- dz k --- dzttruA^-'-A^Zn), (A.19) 

where here and in all formulae below symbol dz means integration over the contour 
C2 such that the points a, z, b on C2 are placed in the counterclockwise order. Using 
the cyclic invariance of the trace we can rewrite this term as the l/n-th of all cyclic 



v * / 



dz dw . 
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z 4 



Figure 4: A typical ladder diagram for the connected correlator of two 
circular Wilson loops on S 2 . In this example, the points Z\, z±, zj belong 
to T(2, 1) and W5, w 2 , are the corresponding points in I(T(2, 1)). The 
points Z2, z 3 , 25, zq belong to T(2, 2) U/(T(2, 2)), and finally w±, W3 are in 
T(1,1)U J(T(1,1)). 



permutations of the points Zi, . . . , z n on the contour, and we get 
w r (C2) = da 1 - / d*i / dz n --- dz k --- dz 2 ti R {A z (z x ) ■ ■ ■ A z (z n )) 



?? 



''2 







rzk+i 


— / d^i 


/ d^n • • • 








J Z\ 


„ e 27Ti 








■ J zidz k ■ 


■1: 









z\ 



dz x I z x dz n --- I zidzk--- I zidz 2 tT R {A z {zi)A z {z 1 z 2 )---A z {z 1 z n )) 



(A.20) 

where in the last line we have changed integration variables Zk = Z\Zk for k = 2 . . . n. 

Now, using (1A.20I) and (1A.15|)(1A.16|)(1A.17I) . we can explicitly evaluate any given 
Feynman diagram D, which defines for us the sets T(2, 2), T(2, 1), T(l, 1), I(k), J(k) 
introduced above. For each point Zf~ G C 2 for A; > 1 we substitute our change of 
variables Zk = z\Zk- Without loss of generality, let us suppose that the point z\ belongs 
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to the set T(2, 1)0. Then for the given ladder diagram D we get 



92d\ 
d R d R > \ 7T J 



D 

2 \ (n+n')/2 



T22 



dw„i 



ri 



dw. 



n'-l 



W2 ^ />e 2lrI r2 re 9 ™" 



T21 



Tn 



dwr 



''2 



X 



Zi - rf/rfw J( i) 



n 



. fcer(2,i) 



Z\Z k - rj/rfw I{k) 

Wi( k )ZiZ k (ziZ k - W!( k )) 



Z\ dz r 
\ 



1 + rf) 2 

Zk + l 

Z\dz k 



X 



/ 



n 

fceT(2,2) 



Z-LZkZlZlik) 



Zxdz 2 x 



n — 



x 



X 



fceT(2,l) fc€T(2,2) fceT(l,l) 

k^l 



(A.21) 



Now is the key step of the computation. First we evaluate the contour integral over 
z\. The integrand in (IA.21I) is a rational function with respect to Zi, so we can evaluate 
the integral by residues. Using that < \z±\ and that \z k \ = 1, we can see that there 
are no residues outside the integration contour \zi\ = r 2 , except the residue at Z\ — 00. 
So, taking the residue at Z\ = 00 we get 



D 



d R d R > \ 7r J 



2 \ (n+n')/2 



— Tn 



T22 



-rf 1 



dw,„> 



n 



n 



"W2 J 

dw n /_i I dwi — 

n ./! 



1 + rf 1 + r\ 

1 



T21 



-r 



Tn 



dz k ■ • ■ / dz2 x 



X 



X 



X 



X 



Z 2 Z 3 ...Z n WiW 2 ■■■w n > 

himi) n 8 ik j m y\_ s ikii( k ) n 5 ikjj(k) T h ■ ■ ■ ^ tr R' T n • ■ ■ T j n > ■ 

fcGT(2,l) fc€T(2,2) fceT(l,l) 



(A.22) 



The remaining integrations are now elementary. The integral over z 2 ,...,z n gives a 



factor ^fa-i)! anc ^ ^he integral over u>i, . . . , w n > gives a factor ^77—. Hence, we finally 



(2-Kif 



(n-1) 



17 Since we are computing the connected correlator, T(2, 1) is non-empty and we can always use 
cyclic invariance so that z\ S T(2, 1). In any case, it will be clear from the computation that it does 
not matter to which set z% belongs. 
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get 



Wg\C 2 )W£\G 



(n') t 



D 



■{2m 



) n+n' \ f 9 2 d \ 

n\n'\ \ 7r / 



2 \ {n+n')/2 



T22 



-r\ 1 



d R d Rl x '"' n\n'\ \ vr J V( 1 + r l)V V 1 + r i 1 + r l 7 V(! + r ?) 2 

^iji(i) n ^iz(fe) n Kink) n ■■ ■ ^ ^ ^ ■ ^ ■ ( a - 23 ) 



X 



feeT(2,i) 



fceT(2,2) 



fcer(i,i) 



This ex pre ssion agrees with the corresponding Feynman diagram in the two-matrix 
model 29([30|. To see the equivalence one needs expressions for the areas Ai,A 2 on 
S 2 written in terms of r 1; r 2 . We denote by A\ the area of the disk on S 2 inside C\ 
(|z| < ri) and by A 2 the area outside C 2 (\z\ > r 2 ). Using z = e^tan | we have 



1 + r\ 
Then 



Ai 
A 



l + rf 



A- At 
A 



1 +r 2 



A- A 2 
A~ 



1 + r 2 



A 



(A.24) 



W R n) {C 2 )W R V{C x ) 



D 



1 l -Mt^ ((A ~ A ' l)M 



T22 



A^V 21 f(A-A 1 )A 1 \ 111 



dRd^nln'l^"^" 1 ' V ^ J \ A 2 J \ A 2 

n n n ^/w^^i---^^^---^- ( a - 25 ) 

fceT(2,l) fc€T(2,2) fceT(l,l) 

Thisprecisely agrees with the corresponding Feynman diagram in the two-matrix 
model H Eo 



DXDY exp 



1 ( X a X a (X a — Y a ){X a — Y a ) Y a Y a 



2g 2 d V ii ' A — Ax — A 2 A 2 



3 — t r fl' e — tr/j e , 
(A.26) 



which has been derived from the zero-instanton sector of 2d YM. 
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